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Abstract 

The decays K~^ —?• and Kl —)• being the theoretically cleanest rare decays 

of mesons, are very sensitive probes of New Physics (NP). In view of the excellent 
prospects of reaching the Standard Model (SM) sensitivity for —)• by the NA62 

experiment at GERN and for —)• by the KOTO experiment at J-PARC, we study 

them in the simplest extensions of the SM in which stringent correlations between these 
two decays and other flavour observables are present. We first consider simple models 
with tree-level Z and Z' contributions in which either Minimal Flavour Violation (MFV) 
or a U{2)^ symmetry is imposed on the quark flavour-violating couplings. We then 
compare the resulting correlations with those present in generic models in which the latter 
couplings are arbitrary, subject to the constraints from AF = 2 processes, electroweak 
and collider data. Of particular interest are the correlations with e'/e and Al —)• 
which limit the size of NP contributions to —>■ tt+i/P and Al —>■ depending on 

the Dirac structure of couplings and the relevant operators. But in MFV models also 
the constraint from Bg —?■ turns out to be important. We take into account the 

recent results from lattice QCD and large N approach that indicate e'/e in the SM to be 
significantly below the data. While in many models the enhancement of e'/e implies the 
suppression of Kl —)• we present two models in which e' je and Kl —)• can be 

simultaneously enhanced relative to SM predictions. A correlation between A"*" —)• 
and B —)• A(A*)^+/i“, found by us in the simple models considered here, should be of 
interest for NA62 and LHGb experimentalists at GERN in the coming years. The one 
with B —)• A(A*)zzp will be tested at Belle II. 
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1 Introduction 


After more than twenty years of waiting [^, the prospects of measuring the branching 
ratios for two golden modes —)• vr+i/p and —?■ within this decade are very 

good. Indeed, the NA62 experiment at GERN is expected to measure the K~^ —)• 
branching ratio with the precision of ±10% [^|^, and the KOTO experiment at J-PARG 
should make a significant progress in measuring the branching ratio for Kl —)• j4j[^. 

These decays are theoretically very clean and their branching ratios have been cal¬ 
culated within the SM including NNLO QGD corrections [6||^ and NLO electroweak 
corrections [9-11 . Moreover, extensive calculations of isospin breaking effects and non- 


perturbative effects have been done 12,13 . Therefore, once the CKM parameters |14fe|, 
\Vub\ and 7 will be precisely determined in tree-level decays, these two decays will offer 
an excellent probe of the physics beyond the SM. Reviews of these two decays can be 
found in [4,14 -17. 
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In a recent paper [18| we have reviewed the status of these decays within the SM 
taking into account all presently available information from other observables and lattice 
QCD. In calculating the branching ratios for these decays we followed two strategies: 

Strategy A: in which the CKM matrix is determined using tree-level measurements 
of 

IK.I, \Vub\, 7, (1) 

where 7 is one of the angles of the unitarity triangle. As new physics (NP) seems to be by 
now well separated from the electroweak scale, this determination is likely not polluted by 
NP contributions allowing the determination of true values of all elements of the CKM 
matrix. Inserting these values into the known expressions for the relevant branching 
ratios allowed us to determine the SM values for these branching ratios independently 
of whether NP is present at short distance scales or not. We found 


B{K~^ — >■ 

= (8.4 ± 1.0) X 10“1\ 

(2) 

B{Kl —> TT^Uu) 

= (3.4 ±0.6) X 10“^^ 

(3) 


This strategy is clearly optimal as it allows to predict true SM values of these branching 
ratios. 

Strategy B: in which it is assumed that the SM is the whole story and the values 
of CKM parameters are extracted from AF = 2 observables, in particular ek, AMg, 
AMd and mixing induced CP asymmetries S^Ks a-iid Having more constraints, 

more accurate values of \Vcb\-, |I4b| and 7 than in strategy A could be found implying 
significantly more accurate predictions 


B{K~^ — >■ 7T~^l/u) 

= (9.1 ±0.7) X 10"^\ 

(4) 

B{Kl —)• -K^vv) 

= (3.0 ±0.3) X 10“^^ 

(5) 


These latter results are useful in the sense that in the case of future measurements of 
these two branching ratios differing from them would signal the presence of NP but this 
NP would not necessarily be contributing to these two decays as it could also pollute the 
determination of CKM parameters through loop decays. 

Evidently, strategy A is superior to strategy B in the context of NP analyses, since 
it allows to determine the CKM matrix elements independently of NP effects which 
may depend on a large number of parameters. But in a given NP model, in which 
contributions to rare processes involve only a small number of new parameters in addition 
to the SM ones, strategy B could also be efficiently used. However, in the present paper 
we will exclusively use the strategy A. 

The decays K~^ —)• and —?■ have been studied over many years in 

various concrete extensions of the SM. A review of the analyses performed until August 

More recent reviews can be found in IMth 


2007 can be found in 14 


19,20 . Most 


extensive analyses have been performed in supersymmetric models MM , the Littlest 
Higgs (LH) model without T-parity 1^, the LH model with T-parity (LHT) [^, 28 


Randall-Sundrum models pOpO] , models with partial compositeness and 331 models 
32-34 . All these models contain several new parameters related to couplings and masses 
of new fermions, gauge bosons and scalars and the analysis of —)• and Kl —)• 

TT^vv requires the inclusion of all constraints on couplings and masses of these particles 
and consequently is rather involved. Moreover, the larger number of parameters present 



















1 Introduction 


3 


in these models does not presently allow for clear cut conclusions beyond rough bounds 
on the size of NP contributions to —)• and Kl —)■ 

Therefore, we think that presently in order to get a better insight into the structure 
of the possible impact of NP on K~^ —)• 7r'*“nn and Kl —)• decays, and in particular 

on the correlation between them and other observables, it is useful to consider models 
with a only small number of parameters. With this idea in mind we will consider: 

• General classes of models based on a U{3)^ flavour symmetry (MFV), illustrating 
them by means of two specific models in which quark flavour violating couplings of 
Z and of a heavy Z' are consistent with this symmetry. 

• Models in which the flavour symmetry U{3)^ is reduced to U{2)^, illustrating the 
results again by means of two simple Z and Z' models. 

• The Z and Z' models with tree-level FCNCs in which the quark couplings are 
arbitrary subject to available constraints from other decays. In particular in this 
case we will include right-handed currents which are absent in MFV and strongly 
suppressed in the simplest U{2)^ models. 


Note that in each case we consider as benchmarks Z and Z' models with tree-level FC¬ 
NCs to quarks, and flavour-conserving, as well as flavour universal, couplings to leptons. 
Neglecting the tiny neutrino masses, one can assume NP to have only left-handed vector 
couplings to the neutrino pair, and ignore scalar currents. Therefore simplified mod¬ 
els involving gauge-bosons form a good generalisation of the more specific NP models 
available. The simplified Z can mimic modified Z penguins for instance, occurring in 
supersymmetric models for example, while a Z'-\\ks, particle occurs in several of the other 
models listed earlier. 

In addition to K~^ —?■ -k^vv and —>• the ratio e'/e belongs to the most promi¬ 

nent observables in iF-meson physics. It is also very sensitive to NP contributions, but is 
unfortunately subject to large hadronic uncertainties present in the matrix elements of 
QCD and electroweak penguin operators. Moreover, strong cancellations between these 
two contributions make precise predictions for e'/e in the SM and its various extensions 
Reviews of e'/e can be found in 


difficult. 


35-39 


The most recent analyses of e'/e 
within Z{Z') and 331 models have been presented in and 1^, respectively. See also 


our SM analysis in 18 


Most importantly, improved anatomy of s'je within the SM have been presented 
. It was triggered by the first result on e'/e from the RCB-UKQCD lattice col¬ 
laboration [42], which indicated that e'/e in the SM could be significantly below the 


m 


data, but the large theoretical uncertainties in this calculation did not yet allow for 
firm conclusions. These uncertainties have been significantly reduced in [41] through the 
extraction of a number of hadronic matrix elements of contributing operators from the 
CP-conserving K —)• yttf data. Parallel to this study an important upper bound for the 
contribution of QCD penguins to e' je has been derived from the large N approach |43j . 
The analysis in [^ combined with the bound in 43 demonstrates that indeed e' je in 


the SM could turn out to be significantly lower than its experimental value. We will be 
more explicit about this in section]^ 

Now, in most extensions of the SM found in the literature the enhancement of e'/e 
through NP usually implies the suppression of the branching ratio for RTl —)• But, 

as we will demonstrate in Section 4.6 simplified models can be constructed in which e'/e 
and the branching ratio for —)• can be simultaneously enhanced over their SM 

values. 
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Our paper is organised as follows. In section]^ we collect basic formulae for —)• 

7r“’“nn and Kl —)• valid in any extension of the SM and discuss their general prop¬ 

erties. In section we formulate the simple Z and Z' models in question. In section 
we recall some aspects of e' je concentrating on the simplified models of the previous 
section. In particular we present two simplified models in which e' je^ B{K~^ —)• vr+nn) 
and B{Ki —?• tt^uv) can be enhanced simultaneously over their SM values. In section]^ 
we present formulae for various decays and observables in the simplihed models of sec¬ 
tion!^ and discuss their correlations with —)• Tr+nn and Kl —)• This includes 

b —>■ transitions, B —)• K{K*)h'u and Kl —)• ■ Kl —)• plays an impor¬ 
tant role in constraining the allowed values of B{K^ —)• While some numerical 

results will be shown already in previous sections the main numerical analysis of the 
models of section is presented in section We conclude in section 


2 General Formulae and Properties 


2.1 General Expressions 


The branching ratios for K~^ —)• 7r'*'nn and Kl —)• in any extension of the SM in 

which light neutrinos couple only to left-handed currents are given as follows 


B{K'^ —^ TX'^Ul/) = K+(l -I- Aem) • 

> 0 /^ImXefry 

B{Kl 7X UU) = kl ■ I - ^^5 - I 


Im Xes V , Ac , . Re Xes \ ^ 


where [13] 


K+ = (5.173 ±0.025) • 10“^^ 
kl = (2.231 ±0.013) • 10“^° 


r A 

0.225 

r A 

0.225 


Aem = —0.003, 


( 6 ) 

(7) 

( 8 ) 
(9) 


and Aj = V*gVid are the CKM factors. For the charm contribution, represented by Pc{X), 
the calculations in [7||^ [T^[T^ imply 

PciX) = 0.404 ± 0.024, (10) 


where the error is dominated by the long distance uncertainty estimated in j^. In 
what follows we will assume that NP does not modify this value, which turns out to be 
true in all extensions of the SM we know about. Such contributions can be in any case 
absorbed into the function X^g. The latter function that describes pure short distance 
contributions from top quark exchanges and NP is given by 

= VZVtd{XL{K) ± Xr{K)) ^ VZVtdXl^{K){l ± (11) 


The functions Xl{K) and XLi{K) summarise the contributions from left-handed and 
right-handed quark currents, respectively. In the SM only Xl{K) is non-vanishing and 
is given by [Is] 


Xf^{K) = 1.481 ± 0.005th ± O.OOSexp = 1.481 ± 0.009. 


(12) 















2 General Formulae and Properties 


5 


One can also express the function Xefj as a function of the branching ratios B{K'^ —)• 
7r“’“nn) and B{Kl —?• vr^nn), which is useful for the study of correlations of the latter with 
other flavour observables. One has, directly from Hi, 0, 




ImXeff 




B{K+ 


TT ' nn 


B{Kl 


0 -mV2 


K+(l + Aem) 


kl 


-X^ReXcPciX), 


B{Kl —)■ TT^VU) 


KL 


1/2 


(13) 

(14) 


In choosing the signs in these formulae we assumed that NP contributions do not reverse 
the sign of SM functions. For more general expressions admitting such a possibility 
see 1441. At the Grossmann-Nir bound (45] the square root in (13) vanishes. 


2.2 Basic Properties 


The correlation between B{K~^ —?■ Tr+nn) and B{Kl —?■ Tr^nn) depends on the short 
distance dynamics, encapsulated in the two real parameters ^ and 9 that vanish in the 
SM. Measuring these branching ratios one day will allow to determine those parameters 
and, comparing them with their expectations in concrete models, obtain insight into 
the flavour structure of the NP contributions. Those can be dominated by left-handed 
currents, by right-handed currents, or by both with similar magnitudes and phases. In 
general one can distinguish between three classes of models 


46 


1. Models with a CKM-like structure of flavour interactions. If based on flavour 


symmetries only, they include MFV and U{2)^ models 47 . In this case the function 


Xl[K) is real and Xji{K) = 0. There is then only one variable to our disposal, 
the value of Xl(K), and the only allowed values of both branching ratios are 
on the green branches in figure But due to stringent correlations with other 
observables present in this class of models, only certain ranges for B{K'^ —)• vr+z/p) 
and B{Ki —?• -k^vv) are still allowed, which we will determine in the context of our 
analysis. 


2. Models with new flavour and CP-violating interactions in which either left-handed 
currents or right-handed currents fully dominate, implying that left-right operator 
contributions to £k can be neglected. In this case there is a strong correlation 
between NP contributions to £k and K —)• ttz^z/ and the £k constraint implies the 
blue branch structure shown in figure On the horizontal branch NP contribution 
to A —)• -Kiyv is real and therefore vanishes in the case of —)• On the 

second branch NP contribution is purely imaginary and this branch is parallel to 
the Grossman-Nir (GN) bound (^. In practice, due to uncertainties in ek, there 
are moderate deviations from this structure which is characteristic for the LHT 
model 27 , or Z or Z' FGNG scenarios with either pure LH or RH couplings |48||49 


3. If left-right operators have signihcant contribution to £k or generally if the cor¬ 
relation between ek and K —>■ ttvv is weak or absent, the two branch structure 
is also absent. Dependent on the values of ^ or 0, any value of B{K^ —>■ vr+z/p) 
and B{Kl —)• is in principle possible. The red region in figure shows the 

resulting structure for a fixed value of ^ and 0 < 9 < 27r. Randall-Sundrum models 
with custodial protection (RSc) belong to this class of models (^. However, it 
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B{K+ ^ TT+uu) 

Figure 1: Illustrations of common correlations in the B{K^ — )• vr+i/p) versus B{Kl — )■ vr^i/p) 
plane. The expanding red region illustrates the lack of correlation for models with general 
LH and RH NP couplings. The green region shows the correlation present in models obeying 
CMFV. The blue region shows the correlation induced by the constraint from Ek if only LH or 
RH couplings are present. 

should be kept in mind that usually the removal of the correlation with ek requires 
subtle cancellations between different contributions to ek and consequently some 
tuning of parameters (^[^ . 

Unfortunately, on the basis of only these two branching ratios alone, it is not possible 
to find out how important the contributions of right-handed currents are, as their effects 
are hidden in a single function In this sense the decays governed by 6 —)> sni? 

transitions, which will also enter our analysis, are complementary, and the correlation 
between K —?■ tti^P decays and B —?■ K{K*)uP, as well as Bg^d j can help in 

identifying the presence or absence of right-handed currents. 


3 Simplified Models 


In studying correlations between various decays it is important to remember that 

• Correlations between decays of different mesons test the flavour structure of cou¬ 
plings or generally flavour symmetries. 

• Correlations between decays of a given meson test the Dirac structure of couplings. 


We will look at the first correlations by comparing those within MFV models based 
on a U{3)^ flavour symmetry with the ones present in models with a minimally broken 
U{2)^ flavour symmetry 50,51 . In the latter case we will work at leading order in 


the breaking of the symmetry, and therefore assume that only the left-handed quark 
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couplings are relevant, as in MFV. We will then extend the analysis to more general 
models with generic flavour structure. 


3.1 Z models with flavour symmetries 

In order to exhibit correlations of —)• and Kl —)• decays with other 

observables we will first consider two simple Z models in which the quark flavour violating 
couplings are consistent either with a U{3)^ or with a U{2)^ symmetry. These models 
are very restrictive as the Z mass and its couplings to leptons are known. In particular, 
in the conventions of |48| for the couplings A(Z) of the Z boson to fermions, 

Af (Z) = A>f{Z) = — = 0.372 . (15) 

2cw 

However, in order to be able to generalise our analysis straightforwardly to the Z' case, 
we will use the general expressions for these lepton couplings. 

We will then find that in the case of MFV there is only one new real parameter a 
and in the U{2)^ case there are three new real parameters: real a and a complex b. 


3.1.1 U{3)^ case 

In this case the Z quark flavour violating couplings are given respectively for the three 
meson systems {K, Ba, Bg) as follows: 

Af (Z) = aVZVtd, Af (Z) = aV:,Vtk, Af (Z) = aV^Ytt, (16) 


where a is flavour-universal and real. 

The presence of tree-level Z contributions in various flavour observables can be sum¬ 
marised by shifts in the master functions S, X and Y which enter respectively the 
expressions for quark mixing (AT = 2) and branching ratios for meson decays with vi? 
and /r"*"//” in the final state. 

The couplings in (16) imply then: 


AS{K) = AS{Bd) = AS{Bs) = AS = 

where 

siM = 4!^ = 1.78137 X 10-' GeV-'^ . 

with Gp being the Fermi constant, f = 1.068 is a QCD correction [48| . 
Similarly, 


4f 

^z9sm 


(17) 


(18) 


Z) 

AXl{K) = AXpiBd) = AXl{Bs) = AX = (19) 


and 


AYa{K) = AYA{Bd) 


AYa{Bs) = AY = a 


K^jZ) ^ 

^z9sm 


( 20 ) 


We observe very strong correlations between the three meson systems. This model has 
only one new real parameter a with respect to the SM, which could be positive or 
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negative. In fact, using the equality of the Z couplings in (15) and eliminating the 
parameter a we find a very stringent relation 


AX = AF = ±4.67\/A5, 


( 21 ) 


where the sign corresponds to two possible signs of a. The consequences of this relation 
are rather profound. In particular: 

• The size of possible effects in rare decays is strongly bounded by the allowed uni¬ 
versal shift in the box function S. 


However, as S'sm > AisM > Fsm > 0, NP generically affects, in this scenario, rare 
decays stronger than particle-antiparticle mixing. 

While the flavour universal shifts AX and AY can have generally both signs, with 
the real parameter a, the universal shifts AS" are strictly positive in agreement with 
the general discussion in 52 . This means that AMg^d and ex can only be enhanced 


in this scenario, and this happens in a correlated manner. 

Due to the present data on Bg —)• the shift AY < 0 is favoured, implying 

suppression of all rare decay branching ratios governed by the functions X and Y. 
Moreover, the amounts of these suppressions are correlated with each other. We 
stress that this property is characteristic for tree-level Z exchange and originates 


in the signs of the leptonic couplings in (15). 


As in the SM Xsm > Tsm, NP affects stronger decays with pp in the final state 
than those with uu. 


Our numerical analysis in section will show that in this scenario NP effects are generally 
below 50% at the level of the branching ratios. 


3.1.2 U{2)^ case 


The Z couplings in (16) are now modified to 


Af{Z) = aVZVtd, Af{Z) = bV^dVth, Af (X) = bVZVtf,, (22) 

with b ^ a being a complex number. Therefore, compared with the U{3)^ case, b 
represents two new real parameters: its absolute value, and the phase which has impact 
on CP violation in Bg^d systems. In this case the correlation between the K system and 
the Bg^d systems is broken. For the K system the MFV formulae remain unchanged, 
while now 


We find then 


A /Y' 

AS{Bd) = AS{Bg) = AS{B) = (6*)^-^, 

(23) 

AXUBd) = AX^Bg) = AX{B) = 

Mz9sm 

(24) 

AYA{Bd) = AYA{Bg) = AY{B) = . 

Mz9sm 

(25) 

AX{B) = AY{B) = ±YQ7^JAS{B)*. 

(26) 
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Moreover, writing the total S function as 


S{B) = SsM + ^S{B) = \S{B)\e-^^^ , 


(27) 


where a non-zero (p is generated by quark flavour violating Z couplings, we find the 
known anti-correlation between mixing induced CP asymmetries in B^ and Bg systems 
respectively: 

S^Ks = sin(2/3 -k 2ip), = sin( 2 |d<i| - 2ip) (28) 

We note then: 


While ek can only be enhanced in this scenario, the fact that 6 is a complex number 
implies the possibility of |5'(-Bq)| being larger or smaller than S'sm, and therefore 
allows for both enhancements and suppressions of AMg^d, independently of ek- In 


this manner some tensions in the unitarity triangle fits can be avoided 150,51 


In the Bg d meson systems the suppressions of branching ratios are favoured by the 
Bg —)• fi~ data. 


Due to the measured value of being SM-like, also the size of allowed modihca- 
tions in is predicted to be small. As seen in (28) the modifications of these two 
asymmetries are anti-correlated with each other and for fixed 7 this anti-correlation 
depends on the value of \Vub\ (^- Similarly to MFV, this scenario favours then 
\Vub\ from exclusive decays, although it still allows for visible non-MFV effects. 


• However, due to the breakdown of the correlation between Bg^d and K meson 
system, NP effects in —>■ vr+z/P and Kl ^ can be larger than in the MFV 

case, being only subject to constraints from ek, AMk, Kl ^ and also e'/e. 

As we will see below the absence of correlation with Bg —?■ is important here. 


3.2 Z' models with flavour symmetries 

These models are less restrictive, and in the MFV case have four new real parameters 
relative to the SM, 

a, Mz', Af(Z'), A7(Z'), (29) 

where A{Z') denote the Z' couplings to fermions, and this number is reduced in the 
correlations between various observables. In the case of U{2)^ models an additional 
complex parameter b ^ a in Bg^d systems is present and the correlations between the K 
system and the Bg^d systems are broken. 

The relevant formulae for the shifts in various functions are obtained from the ones 
in the Z models by simply replacing Mz by Mz' and the Z couplings by Z' ones. The 
QCD correction f in © depends logarithmically on the Z' mass [^. For definiteness 
we will set r = 0.941, which corresponds to Mz' = 5TeV. 

The crucial difference between Z' and Z models is not only the big difference in their 
masses but more importantly that the Z' couplings to leptons are in principle arbitrary 
and do not have to satisfy the relation ©• On the other hand, in accordance with 
the SU{2)l symmetry we have for all Z' models, independently of whether a flavour 
symmetry is imposed, 

Af (Z') = Af{Z'), A(f (Z') = 2Af (Z') + Af{Z'). 


(30) 
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This will have interesting consequences as we will see below. Moreover, these couplings 
and Mz' are constrained by LEP II and present LHC data. 

The correlations between various loop functions in the MFV case have now the struc¬ 
ture 


2\/f Mz'gsu ’ 

(31) 

2-v/f Mz'gsu ’ 

(32) 


and generally, in contrast to (21), /S.X ^ AT. In the U{2)^ scenario these formulae 


apply separately for the loop functions of the K and systems, which generally differ 
from each other. Notice that A5 is always real in the U{3)^ case. 

The following new features relative to the case of Z models should be noted 


As now A'£'{Z') and A^{Z') can differ from each other, the correlations between 
decays with muons and neutrinos in the final state are in principle absent. Therefore 
even in the MFV scenario the data on Bg —?■ alone, being sensitive only to 

A^{Z') have no impact on K~^ —?■ Kl —?• and b —?■ svi' transitions. 

However, when the data on H —)■ fj,~ are taken into account and the 

coupling Ay^(Z') is restricted, the SU{2)i relation in (30) implies some bounds on 
K+ —)■ and Kl ^ in addition to those following from the allowed size 

of AS*. We will be more explicit about this issue in section 5.1 below. 


• After AS* and b —?■ constraints have been imposed, for fixed leptonic Z' 

couplings, NP effects in rare decays decrease with increasing Mz' and as we will 
see in section for Mz' > 5 TeV they will be rather small, in particular smaller 
than in particle-antiparticle mixing. This opposite hierarchy between NP effects 
in mixing and rare decays relative to the Z case could allow one in the future to 
distinguish Z and Z' scenarios. 


• In the U{2)^ scenario also the correlations between NP effects in AT —)• irui' decays 
and Bg d meson systems are broken allowing still significant enhancements of both 
branching ratios subject to the constraints from ek, AMk and the LEP and LHC 
bounds on the Z' mass and its leptonic couplings. 


3.3 Z and Z' with arbitrary FCNC quark couplings 

Finally, we will investigate the cases of general FCNC quark couplings of Z and Z' 
so that non-minimal sources of flavour violation will be present in all meson systems 
and generally they will not be correlated with each other. This will allow larger NP 
contributions to K~^ —)• and Kl —)• than what was possible in the previous 

cases. 

The simplest scenario of NP with non-minimal sources of flavour violation is the 
case of the Z boson with FCNCs. The only freedom in the kaon system in this NP 
scenario are the complex couplings as the Z mass and its couplings to leptons 

are known. In Z' models, in addition to Af^j^{Z'), two new real parameters enter; Mz' 
and A^'^(Z'). In the latter case we will be guided by the bounds on the Z' mass and 
its leptonic couplings from LEP H and the LHC as well as LHCb data on 6 —)• 
transitions. 
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These scenarios have already been considered in [40|48[(54) but the treatment of CKM 
parameters was different there, and both the input from lattice QCD and the value of 
\Vcb\ have changed in the meantime. 


4 e'/e 

4.1 General Structure 

Let us begin our presentation of e'/e with the general formula for the effective Hamilto¬ 
nian relevant for K —)• tttt decays in any extension of the SM 


TiesiK tttt) = HesiK 7 r7r)(SM) + HesiK 7 r7r)(NP) (33) 


where the SM part is given by 

10 

nes{K ^ 7r7r)(SM) = ^ Cf^{lj)Qi (34) 

i=l 

and the NP part by 

10 

nes{K ^ 7r7r)(NP) = ^(G(/x)g, + C'(;u)g'). (35) 

i=l 

Explicit expressions for the operators Qi can be found in [^. For our discussion it will 
be sufficient to have expressions only for the dominant QCD-penguin and electroweak 
penguin operators: 


QCD Penguins: 


Qb = {sd)v-A ^ {mW+A, Qe = {Sadi3)v-A ^ {qi3qa)v+A, (36) 

q=u,d,s,c,b q=u^d^s,Cyb 

Electroweak Penguins: 

^7 = 2 {^d)v-A ^ eg {qq)v+A, Qs = - {sadp)v-A ^ eq{q^qa)v+A- (37) 

q=u,d,s,c,b q=u,d,s,c,b 

Here, a, 13 denote colours and eg denotes the electric quark charges reflecting the elec¬ 
troweak origin of Qj, ..., Qio. Finally, {qq')v±A = 9a7/^(1 =t l^)q'a- The so-called primed 
operators Q[ are obtained from Qi by interchanging V — A and V + A\ these new oper¬ 
ators contribute in the presence of right-handed flavour-violating couplings. Note that 
if NP scales are well above mt, as is the case of Z' models, the summation over flavours 


in (36) and (37) has to include also the top quark. But in the SM and Z models the top 


quark is already integrated out. 

The Wilson coefficients are known at the NLO level in the renormalisation 


group improved perturbation theory including both QCD and QED corrections |55,56 
Also some elements of NNLO corrections can be found in the literature [6,57 


If new operators beyond those present in the SM contribute to e'/e one should in 
principle perform the full RG analysis at the NLO level including these operators. How¬ 
ever, in view of various parameters involved we will follow the procedure proposed in 40 
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and consider NP contributions at the LO. Moreover, as demonstrated there, at the end it 
is a good approximation to include in e'/e only the modifications in the contributions of 
the dominant QCD penguin (Qe) and electroweak (Qs) operators and in the contribution 
of the corresponding primed operators. 

Now, relative to the case of —)• Kl —>■ and A.F = 2 processes, 

flavour diagonal quark couplings are involved, and without knowing these couplings the 
correlation between rare K decays and e'je is lost. In the case of Z the diagonal quark 
couplings are known and this implies a correlation between rare K decays and e’/e, as 
first stressed in [^. But the case of Z' is different. For instance it could be that for some 
reason the flavour-diagonal quark couplings to Z' are very strongly suppressed relatively 
to the non-diagonal ones. In this case one would be able to enhance the branching ratios 
for —7- and Ki ^ without violating the e'/e constraint. We stress this 

point as the usual statements about correlation between rare K decays and e' je made in 
the literature apply to concrete models and one cannot exclude that through particular 
choices of flavour-diagonal Z' couplings to quarks this correlation can be broken. In what 
follows we will restrict our discussion to cases for which such correlations are present. 

Finally, although the impact of s' je also depends on the different scenarios for Z 
couplings, as shown in [^, the SM value of ek must be consistent with the data if one 
wants to satisfy simultaneonsly ek and e'/e. The details depend on the value of the 
hadronic matrix element of the QCD pengnin operator Qg) or equivalently on the value 
of the parameter . If ek in the SM differs significantly from the data, NP required 

to fit the data on ek automatically violates the e'/e constraint for within 20% 

from its large N value = 1.0. But, as we shall see in detail in section 


4.6 


a new 
41 


43 


insight in the range of values of has been gained through the studies in 

so that now more space is left for NP contributions to e'/e. Also, as already mentioned, 
significant arbitrariness in the diagonal quark couplings to Z' allows for larger NP effects 
in this case. 

In |18| we have updated the analysis of e'/e within the SM and the recent analyses 
of e'/e within Z{Z') and 331 models have been presented in |40| and 34 , respectively. 


However, since then two improved analyses of e'/e in the SM have been presented 41,43 
and we will base our analysis on these two papers. 


4.2 SM Contribution 

The starting point of our presentation is the analytic formula for e'/e within the SM 


36,40 , which has been recently updated in 41 and is given as follows 


SM 


= Im [XtFe'{xt )], 


where 


F£/(xt) = Pq + Px Xo{xt) Py Yo{xt) + Pz Zo{xt) + Pe Eoixt) 


(38) 


(39) 


The first term in (39) is dominated by QCD-penguin contributions, the next three terms 


by electroweak pengnin contributions and the last term is totally negligible. The xt 
dependent functions have been collected in the Appendix A of |18|. 
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The coefficients Pi are given in terms of the non-perturbative parameter^ 


Re = 4'^'^ 


114.54 MeV 


1 2 


ms{mc) + md{mc) \ 


Rs = 


114.54 MeV 
msirric) + md{mc) 


as follows: 


(40) 

(41) 


D (0) I (6) D I (8) 73 

Pi = rl + r\ Re + r] Rs • 

The coefficients and comprise information on the Wilson-coefficient func¬ 

tions of the AS = 1 weak effective Hamiltonian at the NLO. Their numerical values for 
three values of as{Mz) are collected in the Appendix B of 1^. We will next describe how 


the (38) is modified in the presence of NP contributions. The structure of modifications 


depends on NP model considered. 

In our numerical analysis we will use for the quark masses the values of 59 , given in 
table Then at the nominal value fx = rric = 1.3 GeV we have 


ms{mc) = (109.1 ± 2.8) MeV, md{mc) = (5.44 ± 0.19) MeV. 


(42) 


4.3 CMFV and U{2f 

These are the simplest cases as only the shifts in the function X, Y, Z, discussed in 
previous section, have to be made if NP is not far from the electroweak scale. The most 
predictive in this case is Z scenario as in this case the following shifts in the functions 


V, Y and Z entering the analytic formula (38) have to be made 


o 2 

AX = AY = AZ = acw^ , 


(43) 


which equal just the shifts in (19). The reason why the shift is universal in these three 


functions originates in the fact that a Z exchange with flavour violating couplings in one 
vertex and known flavour diagonal couplings modifies just the Z-penguin contribution 
which universally enters V, Y and Z. 

The shift AZ has the largest impact on e'/e, as the coefficient Pz is large and negative. 
For a positive a the enhancement of K~^ —)• and —>■ implies suppression of 

e'/e, while a negative a suppresses these branching ratios and enhances e'/e. In fact, in 
MFV this scenario appears to be favoured by the Bg —)• data. Moreover, it would 

also be favoured by the data on e'/e, if the SM prediction for e'/e will turn out to be 
below its measured value, as presently indicated by the analyses in (4l| - (4^ 

The correlation with B{Kl —)• is made manifest using expression (14) together 


with (43), from which one has 


MFV 


+ h{Px + PY + Pz) 


SM 




, (44) 


while the correlation with B{K^ —)• follows from the fact that the phase of X^s 


in (13) is aligned with the SM. 


If the flavour symmetry is reduced down to U{2)^ the formula in (43) is still valid 
but the correlation with Bg^d meson systems is broken and the constraints on the NP 


^Note that Ri do not contain the factor 1.13 present in 40 . 
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contributions to e'/e are weaker. In particular, independently of Bg —)• the ratio 

e'je can be enhanced or suppressed but its MFV correlation with —>• and 

Kl ^ remains valid. 

The case of Z' is complicated by the fact that the diagonal quark couplings are 
rather arbitrary, and are not constrained by other semileptonic rare decays. The analysis 
of e'/e can therefore not be very specific even if constraints from LEP and LHC are 
taken into account. It should also be emphasized that, depending on the structure of 
diagonal couplings, different operators dominate e'je (even if generally they are Qq, 
Qs or the corresponding primed operators). The good news in Z' scenarios is that 
unless a concrete framework is considered, there is no strict correlation between e'/e and 
K —)■ TTi^P allowing for larger NP effects in these decays than what is possible in the case 
of Z scenarios. 


4.4 Z with general flavour-violating couplings 


It should be emphasized, that this scenario can be realized in many models and in the 
case of the absence of a discovery of new particles at the LHC the flavour violating 
couplings of Z could constitute an important window to short distance scales beyond 

the LHC. _ 

For completeness we recall here the formulae for e'/e derived in 


40 . The details 


including derivations can be found there. Relative to the Z' case, discussed subsequently, 
the RG running in this case is simplified by the fact that the initial conditions for the 
Wilson coefficients have to be evaluated at the electroweak scale as in the SM. We 
consider three scenarios for the quark couplings: only left-handed (LH), only right- 
handed (RH), and left-right symmetric (LRS) (^. In the ALRS scenario of the NP 
contributions to Ki^ ^ and —?• vr+i/P vanish and this case is uninteresting from 

the point of view of the present paper. 


4.4.1 LH scenario 

Here the simplest approach is to make the following shifts in the functions X, Y and Z 


entering the analytic formula (38) 40 


XX = XY = XZ = cw 


Stt^ A""* 


(^) 


ff-' 


A* 


(45) 


This formula gives the generalization of the shifts in (43) to arbitrary LH flavour-violating 
Z couplings to quarks. We have then 


LHS 


SM 




(46) 


Z,L 


where the second term stands for the modification related to (|45|). 

Since t 
and B{Kl 


Since the shifts in the loop functions (45) are universal, the correlation between e' 
is again given by (44). On the other hand, since the phase of the 

is lost in this case. 


X^^ coupling is now arbitrary, the correlation with B{K~^ 


4.4.2 RH scenario 


This case is analyzed in detail in section 7.5 in 40 , where it is demonstrated that by 


far the dominant new contribution to e'je comes from the Qg operator. The relevant 



















4 e'le 


15 


hadronic matrix element {Q '^2 = —{Q &)2 and consequently it is known from lattice 
QCD We refer to [40] for details. 


In this case we have then 


RHS 


- + - 
^/ SM V ^ 


(47) 


Z,R 


with the second term given within an excellent approximation by 40 


= -6.2 X 10'^ 


Z,R 


114MeV 


2 

[ 0 ( 3 / 2 ) j 

- 

0.76 


ImAff(Z). 


(48) 


Note that due to the new lattice results in 


msirric) + md{mc) _ 

( 3 / 2 ) 

the central value of has been 


61 


modified relative to 40 where the older value 0.65 extracted from has been used. 

This result implies that ImA|^(Z) must be at most be 0(10“^) in order for e'je to 
agree with experiment. Then, similarly to the LH case just discussed, NP contribution 
to £k are very small and only for CKM parameters for which ek in the SM agrees well 
with the data this scenario remains viable. 

As far as —)• and —)• are concerned we can use the formulae in 48 


Equivalently, in the case of the RH scenario, one can just make a shift in the function 
X{K): 

g 


AX{K) = 


- 1 

> 

s 

\Af{Z)^ 

-dSM^Z- 

[ Ai J 


Anz) = 


2cw 


(49) 


Expressing ImA^ in terms of B{Kl —)• through (14), one then has 


= -32.6 • Rs 


Z,R 




(50) 


4.4.3 General case 

When both Af^{Z) and Aj^(Z) are present the general formula for e'/e is 


+ 


SM 


/\ L 


+ 


/\ R 


(51) 


with the last two terms representing LH and RH contributions discussed above. This 
formula allows to calculate e'/e for arbitrary Z couplings, in particular for the LRS 
scenario where Af^{Z) = A^^(Z), and for the case presented in section 4.6.1 
The numerical analysis of all these scenarios is presented in section [6] 

4.5 Z' with flavour-violating couplings 

We have already emphasized that in general, in the absence of the knowledge of flavour 
diagonal Z' couplings to quarks, there is no correlation between e'/e and K —)• ttvv 
decays. We will therefore not present a numerical analysis of e'/e in Z' scenarios, except 
for one case in section ItHI 

The analysis in 331 models, where the operator Qg turns out to be most important. 


can be found in . On the other hand, in 40 , where the possible impact of Z' on the 
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A/ = 1/2 rule has been considered, the diagonal couplings could be fixed by requiring 
the maximal contribution of Z' to the Aq{K —)• Trvr) amplitude. In this case the operator 
Qq turned out to be most important. As we will see below, a variant of this model turns 
out to be interesting in view of the recent lattice result on e'/e in 42 and recent analyses 


in 41,43 


4.6 Can e'/e and K — nuu be simultaneously enhanced? 

In most extensions of the SM found in the literature the enhancement of the branching 
ratio for Kl —)• through NP usually implies the suppression of s'/e, and vice versa 
an enhancement of s'/e implies a suppression of Kl —)• -k^vu. We have already mentioned 
this feature in the context of our analysis of Z models with MFV after (43). This is 


related to the fact that there is a strong correlation between the negative electroweak 
penguin contribution to s'/s and the branching ratio for —>• Here we would 

like to present two simplified models in which in fact s'/s and B{Kl —)• can be 

simultaneously enhanced with respect to their SM values. 

This case is of interest in view of the recent result from the RCB-UKQCD lattice 
collaboration which indicates that s'/s in the SM could be significantly below the data. 
Indeed, they find in the SM [4^ 

(52) 


(e7e)sM = (1.4±7.0)-I0-^ 


which is by 2.1a below the experimental world average from the NA48 62 and KTeV 


63,64 collaborations. 


(e7£)exp = (16.6 ± 2.3) • 10 


-4 


(53) 


A recent detai’ 


of H 


( 1 / 2 ) 


from 


42 


ed anatomy of s' js in the SM in also confirms that, with the value 
s'js in the SM is indeed significantly smaller than the experimental 
value. Assuming that the real parts of the K —)• tttt amplitudes are fully governed by 
the SM dynamics and including isospin breaking effects the authors of |4I| find 


(e7e)sM = (1.9±4.5)-I0 


-4 


(54) 


which is by 2.9a below (53). Clearly, the size of this suppression of s' js depends sen¬ 
sitively on the value of the dominant source of uncertainty in the prediction of 

s' js in the SM. But even discarding lattice results, and using the recently derived upper 
bounds on Hg and ' from the large N approach 43 , s' js is found typically by a 


factor of two below the data. Motivated by these finding we looked for models in which 
s' js and Kl ^ could be simultaneously enhanced. 


4.6.1 Simplified Z model 

We consider a model in which Z has both LH and RH couplings, but not equal to each 


other, and not differing only by a sign. As seen in (48), in order to obtain a positive 


contribution to £7^ need Im A|/(Z) < 0. But this alone would suppress the rare decay 
branching ratios. The solution to this problem is the contribution of the Qs operator to 
s' js given in (|45|). While this is not evident from this formula, as shown in 40 , for equal 


LH and RH Z couplings this contribution is by a factor of 3.3 smaller than the one in 
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Modified Z, = -0.5 
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e'/e [10-*] 

Figure 2: 95% C.L. allowed regions for e'fe and —)■ Left: model with flavour- 

ehanging Z boson couplings = — 0.5A|^'^. Center: modified Z, LH scenario A|^ = 0. Right: 
5 TeV Z’ with A|j = 1 and = 0.5. The plots are for 3^ = 1 (blue), Bq = 0.76 (green), 
and Bq = 0.57 (red). The hatched regions are the SM predictions at 2a. The gray band shows 
the experimental result for s'/e. 




(48). On the other hand, the branching ratio for —)• is sensitive to the sum of 

LH and RH couplings. Therefore choosing ImA|^'^(Z) > 0 with 


ImAf (Z)| < ImAf (Z) < 3.3|ImAff(Z)| 


(55) 


one can enhance simultaneously e' je and the branching ratio for Kl ^ In doing 

this, have to be kept sufficiently small in order not to spoil the agreement 

of ReAo in the SM with the data. Moreover, the AMk and £k constraints have to be 
satisfied. 

In the left panel of figure^ we show the correlation between e'Je and —)• 

in the case of Aff{Z) = —2A|^Z), and compare it with the opposite correlation that is 
present in the LH scenario (central panel). The different colours correspond to different 
choices of the parameters and 


d(1/2) 

-°6 

= 1.0, 

r(3/2) 

■°8 

= 1.0 

(blue). 

(56) 

r(1/2) 

-°6 

= 0.76, 

r(3/2) 

■°8 

= 0.76 

(green). 

(57) 

d(1/2) 

-°6 

= 0.57, 

r(3/2) 

■°8 

= 0.76 

(red) . 

(58) 


43 


The first choice is motivated by the upper bound from large N approach 

< 1. The second choice uses the central value for from the RBC-UKQCD 




collaboration 


61 


extracted in 


18 


, and assumes that saturating the 


previous bound. Finally, the third choice uses the central values for 
r(3/2) RBC-UKQCD collaboration, with extracted i 


m 


results in 42 


41 


3oth and 

from the lattice 


As expected, in our simple model the requirement of satisfying the data on e'/e 
automatically implies enhanced values of B{Kl —?• tx^vv), while in the LH model, similar 
to the Littlest Higgs model with T-parity 28 , suppressed B{Kl —)■ tx^vv) is predicted. 
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We do not present the correlation between e'/e and —)> 7r+nn as this also involves 

real parts of the new couplings and is more model dependent. 


4.6.2 Simplified Z' model 

Another example of a model in which B{Kl —?■ and e'/e can be simultaneously 

enhanced has been already considered in . In this model, not the electroweak penguin 
operator Qg, but the QCD penguin operator Qg is affected by NP. A tree-level exchange 
of Z' with left-handed flavour violating quark couplings and flavour universal structure 
of diagonal RH quark couplings generates the Qs operator, and through renormalisation 
group evolution also the Qg operator which at the end dominates the NP contribution 
to e'/e. 

Assuming then that Z' has only LH flavour violating couplings one has 

ImAgP = ImC'g(p)((5g(p))o, (59) 


where 


and 


C^irric) = 1.13 


Af{Z')A‘^^{Z') 

AMl, 


— - 4 


TUsifi) nidifJ-) 


{Fk - F^) B 


( 1 / 2 ) 

6 


(60) 


(61) 


Clearly the size of the NP effects depend on the various couplings of the Z' to quarks 
and leptons. The right panel of figure shows the results for the values 


A^^{Z') = 1, 


Al-{Z') = 0.5, 


(62) 


which satisfy the LHC bounds on flavour-conserving four-fermion interactions, and again 

( 1 / 2 ) ( 3 / 2 ) -- 

for the three choices of the parameters ' and Rg of (56)“(58). 


5 Relations to other Ai^ = 1 processes 

5.1 h —S/X+/X 

It is of interest to see how the decays —>■ and Kl —>■ are correlated 

with b —7- transitions and in particular what are the implications of the B —?■ 

K{K*)pA anomalies for —?■ and Kl —)• in the context of the simplest 

models. 

Let us first note that Z models of any kind cannot explain these anomalies for various 
reasons. In concrete models these anomalies are most easily explained through the shifts 
in the Wilson coefficients Cg and Cio of the operators 


Q 9 = 


Qio = 


(63) 


with 165-68 


^NP 


“'-^10 


-(0.5 ±0.2). (64) 

The solution with NP present only in Cg, with ~ —1, is even favoured, but much 
harder to explain in the context of existing models. We refer to 68 for tables with 
various solutions. 
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Z' (5 TeV) : Constrained MFV 



Figure 3: Allowed ranges for versus B{Kl —)■ vr^i/p) (left panel) and B{K~^ —)■ 

TT+i/p) (right panel) in a simplified 5 TeV Z' model obeying CMFV. The 2 a confidence regions 
shown correspond to constraints from kaon mixing (blue), B mixing (red) and b —?• spApT 
transitions (grey) (from (68l). 


This relation is very badly violated in Z models for which one has 

A7(Z) 


Af{Z) 


= -13.3 


(65) 


in drastic disagreement with (64). The explanation of R —)• K*p~ anomalies would 
then imply very strong suppression of 13{Bg — )• p'^p) relative to the SM which disagrees 
with the data. On the other hand the agreement with the data on B{Bs —)• p~^p) would 
allow only very small value of 

In Z' models we have generally 


sin^ 9wCg = — 


sin^ OwC^Q = — 


Therefore for 


Af{Z') Af^iZ') 

VfgVtb 

(66) 

Af{Z') Af{Z') 

VfgVtb 

(67) 

-Af{Z') 

(68) 


the relation between and in (64) can be satisfied. This is the case of Z' with 
purely V — A couplings both in the quark and lepton sector. 

But the SU (2 )l relation in ( p0| then implies that 

Af{Z') = Af{Z'). (69) 


In turn in the case of MFV, when the first ratio on the r.h.s in (66) and (67) reduces to 
flavour independent a, we have 


AXl{K) = AXLiBd) = AXl{Bs) = ax = - sin" 


NP 


(70) 
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AYl{K) = AYL{Bd) = AYl{Bs) = AY = sin^ ■ (71) 

Therefore, for Z' models with MFV quark couplings, the B —?■ fj,~ anomalies 

imply: 

• Enhancement of the branching ratios B{K^ —)• vr+nn) and B{Kl —?■ relative 

to their SM values; 

• Suppression of the branching ratios B{Bs —)• and B{B(i —)■ lA^~) relative to 

their SM values; 

• Enhancement of the branching ratios B{B ^ K*ui)) and B{B ^ Kvu) relative to 
their SM values as already pointed out in [^ . 

The first of these results does not apply beyond MFV, even in U{2)^ models, but the 
second and third remain true in U (2)^ models. Moreover, for arbitrary Z' quark couplings 
the correlations between B —)■ , B{Bs —)• //“'"//“) and B{B —>■ K*uv) exist due 

to the SU{2)l relation in (30) as already known from other analyses, in particular [69|. 


In the latter case we can compare the region still allowed for 5 TeV Z' shown in the right 
panel of figure with the fit results on Cg from 


68 


In figure we show the regions still allowed in the Cg^ = —versus B{Kl 
and B{K^ —)• Tr~^vi)) planes, in a simplified 5 TeV Z' model obeying CMFV. 

We observe that for < —0.3 one leaves the 2cj range allowed by AMg^d^ a-nd for 
< —0.5 the one allowed by ek and AMk- Thus a massive Z' with MFV couplings 
can lower the tension of the theory with data but cannot fully explain the observed 
anomaly. 


5.2 B K{K*)vv 

There are many reasons for performing an analysis of i? —)• decays in our paper: 

• It is well known that they are strongly correlated with K — )• tiuv decays in models 
with MFV [44] , but also in more complicated models 

• As recently shown in [69] these decays, when measured, could allow to distinguish 

between various explanations of the present anomalies in 5 —)• transitions. 

• It should also be stressed that these decays are of interest on its own as they are 

theoretically cleaner than B —)• and allow good tests of the presence of 

right-handed currents and in general of NP. 


15 


Both decays should be measured at Belle II. The most recent estimate of their branch¬ 
ing ratios within the SM reads (6^: 


B{B+ K+ui?) 


0.0409 


(3.98 ±0.43) X 10"®, 


(72) 




0.0409 


(9.19 ±0.86) X 10“®, 


(73) 


where the errors in the parentheses are fully dominated by form factor uncertainties. We 
expect that when these two branching ratios will be measured, these uncertainties will 
be further decreased and \ Vcb\ will be precisely known so that a very good test of the SM 
will be possible. 
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Figure 4: Allowed ranges for B{K^ —)• vr+z/P) versus B{Bd —)■ K*vv) in a simplified Z model 
(left panel) and a 5TeV Z' model (right panel) obeying CMFV. In the left panel the 2 a 
confidence regions shown correspond to constraints from e'/e (green), Ki^ —)■ (yellow) 

and Bg —)• pApT (magenta), while in the right panel they correspond to constraints from kaon 


mixing (blue), B mixing (red) and b —)■ sp^p transitions (grey) (from ^68l). 


An extensive analysis of these decays model independently and in various extensions 
of the SM has been performed in but only the correlation of —)• vr+nP with 

the b —)• snu in MFV can be found in figure 2 of that paper and we would like to 
extend this discussion. In view of the fact that B —)• AThlj/p decays are correlated with 
B RA p'^p in Z and Z’ models and there are also correlations between B —)• rAiii/ 
and K —ttuP decays in such models, we will find correlations between —)• vr+nP, 

Kl —)• Avv and B —>■ which can be tested by LHCb and NA62 before Belle 

will test the correlations between B —)• rAi/d decays and B —)• rA^+^- analyzed in 
detail in (6^ . 

All formulae necessary for our analysis can be found in 69 and will not be repeated 
here (see in particular section 4.1 of that paper). 

In figure]^ we show the regions allowed at 95% C.L. in the B{R^ —?■ 7r'’'z^P) versus 
B{Bd —t R*vv) plane for a simplified Z and a 5TeV Z' model obeying CMFV. We do 
not show corresponding plots for B{B~^ —)• R^nu) because in CMFV the NP dependence 
is the same as for B{Bd —t R*vv). 


5.3 Kl -a 

Only the so-called short distance (SD) part of a dispersive contribution to —)• p^p~ 

can be reliably calculated. It is given generally as follows (A = 0.2252) 

B{Rl ^ p+p-)sD = 2.01 • 10-9 ' , (74) 


where at NNLO 


pfiY) = 0.115 ±0.017. 


(75) 
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The short distance contributions are described by 


Reff = VtlVtd {Yl{K) - Yr{K)) , 


(76) 


with 


Yl^{K) = vvYoixt), 


r]Y = 0.9982, (77) 


also entering B^ d decays. Notice the minus sign in front of Yr, as opposed 

to Xr in 0, that results from the fact that only the axial part contributes. This 
difference allows to be sensitive to right-handed couplings, which is not possible in the 
case of R —)• -kuv decays. 

In the case of tree-level Z exchange we have 


Yl{K) = Yl^{K) + 


Xf{Z) Xf{Z) 
^tlVtd ’ 


Yr{K) 


AfjZ) Ag^(Z) 

aluMl v,iVtd ’ 


(78) 


with analogous expressions for the Z' case. 

If Y{K) is related to X{K), as in most of the models considered here, one can write 
B{Kr —)■ //■'■//“) in terms of B{K^ —)• n^vi') and B{Kl —)■ in analogy to (44), as 

B{Kl p+p") = 2.01 • 10"^ 


*(bsM T -^sm) H- v~^{Pc{Y) =F Pc{X)) 


A5 ’ A 

B{K~^ —>■ B{Kr —)• 


K+ 


I^L 


(79) 


where the first choice of signs holds whenever only left-handed contributions are present 
- i.e. in MFV, U{2)^, and in the LH scenario for generic couplings ~ while the second 
choice holds for RH couplings. NP contributions to B{Ki —)• vanish in the LRS 

scenario for Z and Z'. 

The extraction of the short distance part from the data is subject to considerable 
uncertainties. The most recent estimate gives [71] 

R(RL^/iV)sD<2.5•10-^ (80) 


to be compared with (0.8 ±0.1) • 10“® in the SM. 

As a preparation for the next section it is useful to recall what is the structure of the 
impact on K~^ —)• and Ki —)• of the constraints from e'/e and Kl —)• 


which have an important interplay |40| : P/s puts constraints only on imaginary parts 
of NP contributions while Kr ^ only on the real ones. As demonstrated already 

48 , the impact of the latter constraint on —?■ and Kr —)• depends 


m 


strongly on the scenario for the Z flavour violating couplings. 


6 Results and comparison of bounds 

6.1 Preliminaries 

The detailed phenomenology in the general case of Z and Z' scenarios, including er, 
AMr and rare decays K~^ —)• vr+z^i/, Kr ^ and Kr ^ has been presented 
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iKfel 

3.88(29) X 10-3 

18 


Fk 

156.1(11) MeV 

59 


40.7(14) X 10-3 

18 


Bk 

0.750(15) 

59 

75 

7 

(73.2/?:3)° 

76 


Fb, 

190.5(42) MeV 

59 

IK/ 

0.2252(9) 

77 


Fb. 

227.7(45) MeV 

59 

\^k\ 

2.228(11) X 10-3 

78 


Fbs \J Bbs 

266(18) MeV 

59 

AMk 

0.5292(9) X 10-2 ps-3 

78 


e 

1.268(63) 


59 

AMd 

0.507(4) ps-3 

77 



0.65(20) 

42 

43 

AM, 

17.761(22) ps-3 

77 


r(3/2) 

^8 

0.76(5) 


61 


1.519(5) ps 

77 


Vcc 

1.87(76) 


79 

BBs 

1.512(7) ps 

77 


Vet 

0.496(47) 


80 

ocs{Mz) 

0.1185(6) 

78 


Vtt 

0.5765(65) 


81 

mciruc) 

1.279(13) GeV 

82 


Vb 

0.55(1) 

81 

83 

m,(2GeV) 

93.8(24) MeV 

59 






mrf(2 GeV) 

4.68(16) MeV 

59 






Mt 

173.34(82) GeV 

84 







Table 1: Values of theoretical and experimental quantities used as input parameters. 


in 


48 and generalized to include V/e in 40 


But MFV has not been considered there 
and it will be of interest to see the allowed size of NP contributions in this case. Earlier 
studies of the upper bounds on NP effects in AF = 2 and AF = 1 processes can be 
found in 72,73 . Here we will concentrate on K~^ —)• and Ki^ —?■ decays but 


will also present some results for other decays. The analyses of rare processes in models 
with an U{2)^ flavour symmetry has been already considered in 51,74 , and in 48 in the 
context of Z and Z' scenarios. But our analysis that uses simple models for couplings 
allows a new insight into these models. 

Also, the present analysis uses a different strategy for the CKM parameters than the 
where various scenarios for these parameters have been considered. In what 


40 


one m 

follows we will use the values of the parameters in Q determined in tree-level decays - 
called “strategy A” in |18| - and we will investigate how large NP effects in K~^ —?■ 
and ATl —)• are still allowed when the constraints from ek, AMk, Kl —)• pVpL~, and 

s' je are taken into account. As already described, the latter constraint will be subject to 
significant non-perturbative uncertainties connected to the parameter In spite of 

this, s'/e already has an important impact on the maximal allowed size of the branching 
ratio, not only for Kl —?■ -k^vv but also for —)• vr+i/P. 

In fact the recent progress on the calculation of e'/e in 41 and [43] , reported already 

40l . In the 


in section makes the impact of this ratio on rare decays larger than in 

^ 3 / 2 ) ( 1 / 2 ) 

following we shall use the lattice value = 0.76(5) from 61 , while for B^ we will 

(1/2) (3/2) 

take an average between the new lattice result 42 and the maximal value Bq ' = Bg ' 
allowed by the large N approach 


43 


In table we summarise the values of the parameters used as inputs in our analysis. 
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Modified Z : Constrained MFV 




0 2 4 6 8 10 12 14 

BR(K+ -)■ -K+vv) [10-“] 


Figure 5: The 95% C.L. allowed ranges for B{Kl —)• 7r°i/p) and B{K~^ —)■ tt+z/P) in a simplified 
Z model (left panel) or a 5TeV Z' model (right panel) obeying CMFV. In the case of the 
smaller U{2fi symmetry, the constraints from B processes can he neglected. Note the difference 
in scale between these plots. 


6.2 CMFV and U{2)^ for Z and Z' models 

In the left panel of figure we show the 2 a allowed ranges from current experimental 
constraints for B{Kl —?• tt^vv) and B{K^ —)• vr^i^P) in a simplified Z model obeying 
CMFV. Similarly, in the right panel of the same figure, we show the allowed ranges 
for a simplified Z' model with a Z' mass of 5TeV, as discussed earlier, also obeying 
CMFV. Neglecting the constraints from the systems gives the situation in the less 
constrained U{2fi symmetry scenario. In both cases we have used the averaged CKM 
inputs from strategy A. We make the following observations: 

• For the simplified Z model, constraints from AF = 1 processes dominate over 
AF = 2 ones. The latter in fact hardly constrain these branching ratios at all. 

• For Z' models the situation is the opposite: due to a direct dependence on the 
high NP scale, AF = 2 observables become the most constraining, and we have 
therefore neglected the AF = 1 constraints. 

• NP contributions in simplified Z models with CMFV are rather constrained by the 

Bg —)■ pFpL~ branching ratio. In Uflfi this constraint is not present, while the short 
distance part of Rl P~ still leaves ample room for NP. On the other hand, 

the strongest limit for an enhancement of B{K^ —)• 7r'*'z^P) and B{Kl —)• Fov) 
branching ratios, both in U{3fi and Uflfi Z models, comes from e'/e. Indeed, 
already the SM point is only marginally compatible with the experimental data, 
and lower values of the two branching ratios are preferred. 

• For Z' models the AF = 2 constraints from the kaon and B systems are comparable 
in size, therefore there is little difference between the CMFV and Uflfi scenarios. 
For a 5 TeV Z' they can deviate from the SM by at most 10 — 20%, which could be 
hard to detect even in the flavour precision era. 

In summary we find that it will not be easy to distinguish MFV models from the 
SM on the basis of —)• vr+i/P and Kl —?• Fuv. While in the case of Z' models 
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b: 

T 

C? 

« 



I LH scenario 
I RH scenario 
I LRS scenario 
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BR(/t+^ ;r+vv) [10~"] 


Modified Z : e'/e, BR(Ki^ AM^; constraints 



Figure 6: The allowed ranges for B{Kl —)■ and —>■ in a simplified Z model 

(left) and a 5TeV Z' model (right) in LH and RH scenarios. The Sk and AM^ constraints 
are imposed in all cases. In the left-handed plot the e'/e and Kl ^ pip. constraints are also 
imposed. 


small NP effects are required by AF = 2 constraints, because of the high Z' mass, in 
the case of Z models the crucial limit comes from the data on Bg —)• p~ and e'/e. 
While an enhancement of the two branching ratios is always strongly constrained, their 
suppression with respect to the SM prediction is still possible in the latter case. 

6.3 Generic Z models 

In the left panel of figurej^we show the 95% C.L. allowed ranges for B{Kl —)• and 

B{K'^ —)• TT^uv) in the LH, RH and LR scenarios with Z mediated FCNC. The origin for 
the different ranges is explained in detail in [40] . Here we only note the following basic 
features: 

• In the LH scenario B{K^ —?• can be by a factor of two larger than its SM 

value. The strong e'/e constraint, on the other hand, forces B{Kl —?■ to be 

of the order of the SM value or smaller, as explained in section [T^ Both branching 
ratios can also be significantly suppressed. We show the impact of the e'/e and 
Kl —)■ p~^p~ constraints. 

• In the RH scenario B{Kl —)• tt^uu) is again constrained to be close to its SM value, 

while B{K^ —)• can be almost by a factor of five larger than its SM value 

because the Kl ^ P^P~ constraint is weaker. Such a large enhancement is anyhow 
already constrained by the present experimental results. Both branching ratios can 
also be suppressed relative to SM values but not as strongly as in the LHS case. 

• Finally in the LRS case the allowed range for B{Kl —>• is similar to the RHS 

case, while, due to the absence of the Kl —)■ p'^p~ constraint, B{K~^ —)■ vr+z^P) can 
be large. The ex constraint plays a role here because of the presence of left-right 
operators. 
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6.4 Generic Z' models 


Due to the sensitivity of the e' je constraint to Z' diagonal quark couplings, in order to 
be model independent, we present a numerical analysis in Z' scenarios without the ^ je 
constraint. In the right panel of figure we show the 3 a allowed ranges for B{Kl —)■ 
and —)• in a simplified 5TeV Z' model for the LH scenario obeying 

the ek and AMk constraints. The leptonic Z' couplings have been fixed to the Z boson 
values for concreteness, = A'£'{Z). Since the AF = 2 effects due to RH currents 

alone are identical to the ones of LH currents, exactly the same results hold also in the 
RH scenario. 

In the LRS scenario the constraints from ex are much stronger, due to the presence 
of left-right operators. Notice, on the other hand, that one can in principle avoid the 
strong AF = 2 bounds by means of some fine-tuning if the RH couplings are sufficiently 


small 49 ; we do not analyse this possibility here. 


7 Summary and Outlook 


In the present paper we have made another look at —>■ and K]^ —)• decays 

which are expected to become the stars of flavour physics in the coming ten years. 
Our results are presented in numerous plots which should allow to monitor efficiently 
the experimental developments in the coming years. In particular the correlations with 
other observables like —)• B —)■ and B —)■ K{K*)ui' branching 

ratios and e'/e will be very relevant for the distinction between various extensions of the 
SM. Also the improvement in the accuracy of the CKM parameters determined in tree- 
level decays and more accurate values of various non-perturbative parameters obtainted 
by lattice QCD will be important ingredients in future analyses. 

In view of the recent result on e' je from RBC-UKQCD collaboration and the 


analyses in [S,43 which find s' je significantly below the data, we have presented two 
simplified models which would improve the agreement of the theory and data if the 
present status of e' je will be confirmed by more precise lattice QCD calculations one 
day. 

We close our paper with the following observations: 

• There is a hierarchy in the size of possible NP effects in 77 —)• 'kvv mediated by 
tree-level Z exchanges. They are smallest in CMFV, larger in U{2)^ models and 
significantly larger in the case of new sources of flavour and CP violation beyond 
these two CKM-like frameworks. 




^ anomalies favour the 
AF = 2 observables however put 


In Z' models with MFV the present B^ 
enhancement of K~^ —)• and 

significant constraints on this possibility. 

Due to the absence of correlation between K —?■ TTvi? and e' je in general Z' models, 
the size of NP contribution in these decays could be large. Then, as demonstrated 

can probe energy scales as large as 1000 TeV 


m 


49 , 77+ ^ 7r+ uu and 


in the presence of general flavour-violating couplings. 

If the NA62 experiment will find the branching ratio for 77+ —)• vr+z^T' to be sig¬ 
nificantly above the SM predictions, both tree-level Z and Z' exchanges could 
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be responsible for these effects - but the same can be said about more compli¬ 
cated models like LHT, RSc and supersymmetric models. Such high values of 
B{K'^ —7- TT^uv) will also signal non-MFV sources at work. 

• In particular, only Z and Z' models with general flavour violating couplings, among 
the models that we considered, allow for B{K^ —)• 'k'^vv) above 20 x 10“^^. 

• Finally, the future measurement of B{Kl —>■ will signihcantly facilitate the 

distinction between various models. 
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